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19.  ABSTRACT  (continued) 


This  approach  provides  a  calculation  scheme  for  studying  guided  wave 
propagation  over  extended  distances,  at  frequencies  of  1  Hz  and  above. 
The  heterogeneity  models  which  have  been  used  are  two-dimensional  and 
calculations  are  carried  out  for  one  frequency  at  a  time. 

A  sequence  of  models  with  varying  levels  of  heterogeneity  have  been 
considered  in  order  to  determine  the  merits  and  limitations  of  the 
computation  scheme.  The  coupled  mode  technique  works  well  with 
heterogeneous  models  in  which  the  local  seismic  velocities  differ  from  the 
stratified  reference  model  by  up  to  2  per  cent  and  there  are  no 
significant  distortions  of  the  main  discontinuities  (e.g.  the  crust-mantle 
boundary) .  The  approach  can  be  used  for  higher  levels  of  heterogeneity  and 
with  distorted  •interfaces  but  a  large  number  of  modes  needs  to  be 
considered  with  consequent  high  computation  co.sts .  If  the  level  of 
heterogeneity  is  not  too  large  then  the  interaction  between  modes  can  be 
restricted,  rather  than  extending  over  the  whole  mode  set,  with  consequent 
reduction  in  computation  co3t . 

One  of  the  major  effects  of  crustal  heterogeneity  is  to  introduce  the 
possibility  of  smearing  out  the  main  amplitude  peak  in  the  Lg  wave  train 
over  a  band  of  group  velocities.  As  a  result,  an  effective  measure  of  the 
energy  content  of  the  Lg  waves  will  be  to  consider  the  integrated 
amplitude  along  the  traces  beween  group  velocities  of  3.6  and  3.3  km/s. 
The  effects  of  heterogeneity  vary  between  different  parts  of  the  Lg  wave 
train  and  the  representation  of  the  wavefield  in  terms  of  modal 
contributions  allows  a  detailed  analyis  in  terms  of  the  group  velocity 
components . 

The  computational  scheme  is  currently  being  extended  to  try  to  allow 
for  mode-coupling  in  three-dimensions  for  both  Love  ana  Rayleigh 
components  of  Lg  waves.  When  heterogeneity  transverse  to  the  propagation 
path  is  relatively  weak  a  good  approximation  to  the  propagation 
characteristics  can  be  produced. 
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SUMMARY 

\ 

'  The  Lg  wave  phase  which  is  of  considerable  interest  for  nuclear  discrimination  problems 
is  normally  observed  after  propagation  through  considerable  distances.  This  phase  is 
dominantly  guided  in  the  crustal  waveguide,  which  is  known  to  be  a  region  with  very 
considerable  horizontal  variability  in  properties. 

The  effect  of  heterogeneous  crustal  structures  on  Lg  waves  has  been  determined  by  using 
a  ’’coupled-mode*  technique  in  which  the  local  seismic  wavefield  in  the  real  medium  is 
expressed  as  a  combination  of  the  modal  eigenfunctions  of  a  stratified  reference  structure. 
Departures  of  the  seismic  properties  in  the  medium  from  those  of  the  reference  medium  lead  to 
coupling  between  the  amplitude  coefficients  in  the  modal  expansion.  The  evolution  of  these 
modal  weighting  factors  with  horizontal  position  are  described  by  a  coupled  set  of  ordinary 
differential  equations. 

This  approach  provides  a  calculation  scheme  for  studying  guided  wave  propagation  over 
extended  distances,  at  frequencies  of  1  Hz  and  above.  The  heterogeneity  models  which  have 
been  used  are  two-dimensional  and  calculations  are  carried  out  for  one  frequency  at  a  time.  — 

A  sequence  of  models  with  varying  levels  of  heterogeneity  have  been  considered  in  order 
to  determine  the  merits  and  limitations  of  the  computation  scheme.  The  coupled  mode 
technique  works  well  with  heterogeneous  models  in  which  the  local  seismic  velocities  differ 
from  the  stratified  reference  model  by  up  to  2  per  cent  and  there  are  no  significant  distortions 
of  the  main  discontinuities  (e.g.  the  crust-mantle  boundary).  The  approach  can  be  used  for 
higher  levels  of  heterogeneity  and  with  distorted  interfaces  but  a  large  number  of  modes  needs 
to  be  considered  with  consequent  high  computation  costs.  If  the  level  of  heterogeneity  is  not 
too  large  then  the  interaction  between  modes  can  be  restricted,  rather  than  extending  over  the 
whole  mode  set,  with  consequent  reduction  in  computation  cost. 

One  of  the  major  effects  of  crustal  heterogeneity  is  to  introduce  the  possibility  of  smearing 
out  the  main  amplitude  peak  in  the  Lg  wave  train  over  a  band  of  group  velocities.  As  a  result, 
an  effective  measure  of  the  energy  content  of  the  Lg  waves  will  be  to  consider  the  integrated 
amplitude  along  the  traces  beween  group  velocities  of  3.6  and  3.3  km/s.  The  effects  of 
heterogeneity  vary  between  different  parts  of  the  Lg  wave  train  and  the  representation  of  the 
wavefield  in  terms  of  modal  contributions  allows  a  detailed  analyis  in  terms  of  the  group 
velocity  components. 

The  computational  scheme  is  currently  being  extended  to  try  to  allow  for  mode-coupling 
in  three-dimensions  for  both  Love  and  Rayleigh  components  of  Lg  waves.  When 
heterogeneity  transverse  to  the  propagation  path  is  relatively  weak  a  good  approximation  to  the 
propagation  characteristics  can  be  produced. 
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RESEARCH  OBJECTIVES 


The  aim  of  this  work  is  to  develop  techniques  for  handling  guided  wave  propagation  in 
three-dimensionally  heterogeneous  media,  especially  for  Lg  waves,  in  order  to  improve  the 
understanding  of  the  nature  of  regional  seismic  phases  and  the  way  in  which  the 
characteristics  of  the  seismic  source  can  be  modified  by  propagation  to  the  receiver.  Such 
information  will  be  valuable  in  assessing  the  effect  of  geological  structure  on  the  behaviour 
of  potential  discriminants  between  earthquakes  and  underground  nuclear  explosions. 

Most  of  the  regional  phases  of  interest  for  nuclear  discrimination  problems  are 
observed  after  propagation  through  considerable  distances.  These  phases  travel  through 
the  crust  and  the  uppermost  mantle,  which  from  a  variety  of  studies  are  known  be 
regions  with  very  considerable  horizontal  variability  in  properties. 

What  is  therefore  needed  is  a  computational  procedure  which  will  allow  the  tracking  of 
guided  seismic  wave  propagation  through  a  horizontally  heterogeneous  crust  for  distances 
up  to  1000  km  or  more.  Such  a  method  should  be  able  to  be  used  with  three-dimensional 
variations  in  seismic  properties  superimposed  on  the  normal  increase  in  seismic  wave 
velocities  with  depth.  A  suitable  candidate  is  the  coupled  mode  scheme  described  by 
Kennett  (1984)  which  has  already  been  used  with  some  success  in  understanding  the 
characteristics  of  Lg  wave  behaviour  (Kennett  &  Mykkeltveit  1984). 


RESEARCH  STATUS  AND  ACHIEVEMENTS 

During  the  past  year,  the  work  on  Lg  phases  has  mostly  been  carried  out  using  the 
coupled  modes  method  introduced  by  Kennett  (1984).  Firstly,  the  limitations  of  the 
calculation  scheme  have  been  explored  for  two-dimensional  heterogeneity  models  and  a 
number  of  approximations  have  been  explored  with  the  aim  of  improving  computational 
performance.  Secondly,  an  attempt  has  been  made  to  extend  the  methods  to  more  general 
three-dimensional  models  of  heterogeneity. 

For  regional  S  wave  trains,  a  representation  in  terms  of  a  limited  number  of  discrete 
modes  gives  a  economical  computational  description  for  horizontally  stratified  media.  At 
each  frequency  only  a  limited  number  of  such  modes  need  be  considered.  In  a  laterally 
heterogeneous  medium,  at  a  single  frequency,  it  is  possible  to  represent  the  seismic 
displacement  and  traction  fields  within  the  varying  medium  as  a  sum  of  contributions  from 
the  modal  eigenfunctions  of  a  reference  structure,  with  coefficients  which  vary  with 
position.  The  evolution  of  these  modal  expansion  coefficient  terms  with  horizontal  position 
can  be  described  by  a  set  of  coupled  partial  differential  equations  in  the  horizontal 
coordinates  (an  outline  of  the  theoretical  development  is  given  for  completeness  in  the 
Appendix).  The  cross-coupling  terms  between  different  modes  depend  on  the  departures  of 
the  heterogeneous  structure  from  the  stratified  reference  model.  These  differences  are  not 
required  to  be  very  small  but  must  not  be  such  as  to  completely  change  the  nature  of  the 
crustal  wave  guide.  Thus,  it  is  possible  to  accommodate  substantial  localised  change  in 
seismic  velocities  and  density,  but  more  difficult  to  allow  for  shifts  of  more  than  2-3  km  in 
the  position  of  major  interfaces  such  as  the  crust-mantle  boundary.  Maupin  &  Kennett 
(1987)  give  an  extended  discussion  of  the  circumstances  in  which  the  coupled  mode 
approach  can  be  applied  to  two-dimensional  heterogeneous  models. 

Studies  of  Two-dimensional  Heterogeneity  in  the  Crust  and  Mantle 

When  the  heterogeneity  within  the  medium  is  two-dimensional,  the  calculations  can  be 
recast  as  the  solution  of  non-linear  differential  equations  of  Ricatti  equation  for  the 
reflection  and  transmission  matrices  connecting  the  modal  expansion  coefficients  at 
different  positions.  The  advantage  of  this  rearrangement  is  that,  for  each  frequency,  the 
boundary  conditions  on  the  differential  equations  are  simplified  for  a  generally 
heterogeneous  medium. 
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If  we  adopt  a  reference  medium  which  does  not  vary  with  horizontal  position,  the 
individual  mode  contributions  propagate  independently  in  that  structure.  However,  once 
the  properties  of  the  true  medium  differ  from  the  reference,  the  independence  of  modal 
propagation  is  lost.  The  effect  of  the  heterogeneity  enters  into  the  differential  equations  for 
the  modal  expansion  coefficients  via  a  coupling  matrix  whose  dimensions  are  dictated  by 
the  number  of  modes  included  in  the  calculation. 

The  choice  of  the  number  of  modes  is  of  considerable  importance.  All  the  significant 
wavenumber  components  for  the  seismic  phase  of  interest  should  be  included,  as  well  as  an 
allowance  for  steeper  angles  of  propagation  than  would  be  present  in  the  reference  medium, 
in  order  to  allow  for  scattering  effects.  The  demands  of  computation  are  that  the  number  of 
modes  should  be  kept  as  small  as  possible,  since  the  computation  time  depends  on  the 
square  of  the  number  of  modes.  On  the  other  hand  if  the  truncation  is  too  tight,  a  poor 
representation  of  the  displacement  field  can  ensue  and  the  results  are  of  limited  use 

In  figure  1  we  illustrate  the  reference  model  (ARANDA)  used  for  most  of  the  Lg  mode 
calculations,  and  the  modal  eigenfunctions  for  the  22  Love  modes  included  in  calculations 
at  1.5  Hz.  The  cut-off  criterion  was  based  on  the  effective  depth  of  penetration  of  the 
modes.  The  fundamental  mode  and  first  15  higher  modes  are  confined  to  the  crust,  and 
are  the  main  contributors  to  Lg.  Mode  16  begins  to  have  significant  displacement  in  the 
mantle,  and  the  remaining  modes  carry  their  main  energy  in  the  mantle  and  represent  the  Sn 
phase.  The  truncation  is  made  at  a  phase  velocity  of  4.54  km/s,  corresponding  in  this 
model  to  a  penetration  depth  of  about  100  km.  Features  in  the  wave  propagation  process 
corresponding  to  higher  phase  velocities  than  4.54  km/s  cannot  be  represented  by  these 
coupled  mode  calculations. 

In  order  to  give  a  good  account  of  the  guided  wave  propagation,  we  must  make  sure 
that  all  the  significant  wavenumber  components  for  the  seismic  phase  of  interest  are 
included,  as  well  as  making  an  additional  allowance  for  steeper  angles  of  propagation  in 
order  to  allow  for  scattering  effects.  This  means  that  we  cannot  just  use  those  modes 
which  have  the  bulk  of  their  energy  within  the  crust  but  must  also  include  waves  with 
higher  phase  velocity  whose  energy  in  the  stratified  reference  model  lies  dominantly  in  the 
upper  mantle.  At  1.5  Hz  as  shown  in  figure  1,  the  fundamental  mode  and  first  15  higher 
modes  are  confined  to  the  crust,  and  are  the  main  contributors  to  Lg.  However,  we  need  to 
carry  at  least  the  22  modes  illustrated  if  we  are  to  allow  for  wave  interaction  processes 
generating  waves  travelling  at  angles  up  to  40’  from  the  horizontal  in  the  midcrust  and  even 
then  we  will  not  give  a  full  account  of  waves  travelling  at  steeper  angles  to  the  horizontal. 

An  indication  of  the  way  in  which  the  coupled  mode  scheme  is  likely  to  behave  when 
confronted  with  heterogeneity  is  provided  by  looking  at  the  expansion  of  the  modal  eigen¬ 
functions  of  a  particular  stratified  structure  in  terms  of  the  eigenfunctions  of  a  different 
reference  structure.  The  situation  is  illustrated  in  figure  2  for  two  cases:  firstly,  where  the 
crust-mantle  interface  is  maintained  fixed  but  the  crustal  velocity  is  varied,  and  secondly 
when  the  crust-mantle  interface  is  moved  by  I  km.  The  pattern  of  modal  behaviour  is 
rather  different  in  the  two  situations.  For  a  change  in  crustal  velocity  (fig  2a),  the  modal 
eiigenfunction  shapes  remain  similar  in  the  varied  model  and  so  the  expansion  coefficents 
are  largest  for  the  original  mode  and  its  immediate  neighbours  and  decay  quite  rapidly  with 
increasing  mode  separation.  A  l  km  shift  in  the  crust-mantle  bounary  has  a  more 
pronounced  effect  (fig  2b),  now  the  modal  expansion  coefficients  alternate  in  sign  and 
decay  less  rapidly  with  increasing  mode  separation.  This  arises  because  the  structure  of  the 
dispersion  equations  and  therefore,  the  nature  of  the  modal  eigenfunctions,  is  much  more 
sensitive  to  the  position  of  a  major  discontinuities  than  to  the  details  of  the  velocity 
distribution  between  them.  A  consequence  is  that  for  a  perturbed  interface  many  modes  are 
needed  in  a  modal  expansion  to  give  an  accurate  representation  of  the  displacement  and 
traction  fields. 


Figure  1:  ARANDA  reference  model  and  Love  mode  eigenfunctions  with  phase  velocity  less 
than  4.54  km/s  at  frequency  1  Hz. 
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Figure  2:  Illustration  of  the  expansion  of  the  eigenfunctions  of  a  varied  model  in  terms  of  those  for  a 
reference  structure.  The  expansion  coeffiecients  are  shown  for  a  span  of  8  modes.  The  reference  model 
and  the  varied  model,  indicated  by  the  chain  dotted  line,  are  shown  in  the  insert, 
a)  change  in  velocity,  b)  change  in  depth  to  crust-mande  boundary 
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Significant  perturbations  of  the  major  interfaces  in  the  model  thus  lead  to  extensive 
coupling  between  modes  and  can  induce  relatively  steep  angles  of  propagation.  If  such 
features  are  required  in  the  models  to  be  considered,  then  a  broad  sweep  of  modes  must  be 
taken  for  the  reference  model  (for  the  Lg  case  at  1.5  Hz,  at  least  30).  With  the  benefit  of 
hindsight,  it  would  appear  that  the  calculations  of  Kci.uett  &  Mykkeltveit  (1984)  on  the 
influence  of  graben  structures  lie  at  the  limit  of  acceptable  model  variations  for  the  modal 
set  adopted. 

The  coupled  mode  representation  is  well  suited  to  heterogeneous  models  where  the 
behaviour  consists  of  relatively  random  variations  about  the  properties  of  the  reference 
model  (see  e.g.  fig  3).  In  this  case,  a  single  incident  mode  will  interact  with  a  limited 
number  of  neighbouring  modes  so  that  the  coupling  matrix  is  diagonal  dominated  with  a 
limited  effective  bandwidth. 

With  the  aid  of  the  coupled  mode  technique,  a  study  has  been  made  of  the  way  in 
which  the  modal  field  for  Lg  waves  is  affected  by  various  levels  of  distributed 
heterogeneity  for  frequencies  up  to  2Hz  over  paths  of  1000  km  long.  A  sequence  of 
different  heterogeneity  models  were  constructed  by  specifying  the  seismic  velocities  in  a 
vertical  section  at  40  km  intervals  and  then  using  bicubic  spline  interpolation  within  each 
layer.  For  the  distributed  heterogeneity  models  the  velocity  values  were  generated  with  a 
random  perturbation  to  the  reference  value  within  a  prescribed  range  of  variation.  The 
velocity  values  were  then  smoothed  horizontally  by  applying  a  moving  average  over  3 
points.  This  imposes  a  typical  horizontal  scale  length  of  around  100  km,  and  the  vertical 
scale  varies  within  the  model  becoming  larger  in  the  upper  mantle.  The  velocity  distribution 
is  arranged  to  return  to  the  reference  structure  at  the  ends  of  the  model  (0  and  1000  km)  so 
that  a  direct  physical  interpretation  of  the  node  coupling  results  in  terms  of  reflection  and 
transmission  can  he  made.  In  figure  3  we  illustrate  a  set  of  models  (D,  E,  F)  with 
increasing  crustal  heterogeneity.  Model  D  is  constricted  to  have  about  ±1  per  cent 
heterogeneity  in  the  crust  with  a  horizontal  scale  length  of  about  100  km,  and  ±1  per  cent 
heterogeneity  has  been  introduced  in  the  mantle  as  well.  Models  E  and  F  retain  the  same 
mantle  heterogeneity  structure  as  model  D  but  have  higher  crustal  heterogeneity  amplitudes, 
±2  per  cent  for  model  E  and  ±5  per  cent  for  model  F.  We  shall  also  consider  model  C  wi 
the  same  crustal  structure  as  D,  but  no  heterogeneity  below  the  crust-mantle  boundary. 

After  propagating  a  large  horizontal  distance  in  a  heterogeneous  model  the  energy 
originally  in  a  single  incident  mode  will  no  longer  be  confined  to  that  mode,  and  indeed 
some  energy  may  be  reflected  back  by  the  heterogeneity,  with  again  the  possibility  of 
conversion  between  modes.  When  the  level  of  heterogeneity  is  relatively  low  (of  the  order 
of  1  per  cent  deviations  from  the  stratified  reference  model  -  as  in  model  D)  transmission 
effects  dominate  and  the  level  of  reflected  energy  is  negligible.  Even  with  rather  larger 
levels  of  heterogeneity  the  reflection  from  the  whole  model  is  small  but,  as  pointed  out  by 
Maupin  &  Kennett  (1987),  it  is  necessary  to  retain  apparently  reflected  waves  in  the  course 
of  the  computation  in  order  to  get  an  accurate  calculation  of  transmission  effects.  The 
Ricatti  equations  are  much  simpler  when  reflection  can  be  neglected  and  so  computation 
time  can  be  reduced. 

For  any  particular  structure,  it  is  therefore  necessary  to  undertake  a  preliminary 
calculation  including  both  reflection  and  transmission  effects  before  it  can  be  assumed  that 
reflection  can  be  ignored.  The  control  of  numerical  accuracy  in  the  course  of  the 
calculations  is  rather  difficult,  although  for  perfectly  elastic  models  there  is  the  possibility 
of  monitoring  the  constancy  of  the  total  energy  in  reflection  and  transmission  associated 
with  an  individual  incident  mode.  To  preserve  full  numerical  precision  quite  small  steps 
have  to  be  taken  in  the  horizontal  direction:  at  1.5  Hz  the  step  should  be  no  larger  than  0.5 
km  if  both  reflected  and  transmitted  waves  are  considered.  For  a  number  of  models  a 
somewhat  larger  step  length  gave  reasonable  results  for  the  case  of  transmission  alone. 
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Figure  4:  Representation  of  the  transmission  matrix  at  I  Hz  after  passage  through  1000  km  of  model  D. 
Amplitude  values  less  than  1  are  left  blank. 
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For  an  incident  mode  at  0  km  the  effect  of  1000  km  of  heterogeneity  is  that  the  energy 
originally  in  the  mode  is  no  longer  confined  to  that  particular  mode.  In  figure  4  we  show 
the  transmission  matrix  for  the  15  mode  set  appropriate  to  the  ARANDA  reference  model  at 
1  Hz,  after  passage  through  the  heterogeneity  model  D.  Each  column  of  the  matrix 
corresponds  to  the  incidence  of  a  single  mode  of  amplitude  100  at  0  km.  We  see  that  the 
behaviour  is  diagonal  dominated,  but  that  significant  mode  conversions  can  occur.  For  the 
Lg  type  modes  (0-9)  the  bandwidth  for  significant  interaction  is  typically  ±  2  modes.  The 
strong  interaction  of  modes  1  and  2  arises  because  of  their  similar  shape  in  the  top  10  km. 
For  the  Sn  modes  (11-15),  the  eigenfunctions  are  very  similar  and  there  is  strong 
interaction  due  to  mantle  perturbation  extending  to  200  km  (the  lower  portion  is  not  shown 
in  fig  3).  However  the  separation  from  the  Lg  modes  is  striking;  there  is  very  little 
interaction  except  for  mode  10  which  shares  some  of  the  characteristics  of  each  group. 

As  the  frequency  increases  the  number  of  modes  need  to  cover  the  same  phase  velocity 
range  increases,  which  causes  some  problems  in  combining  the  results  from  different 
frequencies.  The  increased  size  of  the  differential  equation  system  also  means  an  increase 
in  computation  time  as  the  square  of  the  number  of  modes  if  all  mode  interactions  have  to 
be  considered.  Fortunately,  present  tests  indicate  that  for  the  same  model  the  bandwidth  of 
interaction  increases  only  slowly  with  frequency.  For  the  Lg  waves  propagating  through 
models  similar  to  D  with  about  1  per  cent  heterogeneity,  computations  can  be  restricted  to  a 
band  of  5  modes  either  side  of  the  target  mode  without  appreciable  error.  At  1.5  Hz  with 
22  modes  included,  this  bandwidth  restriction  has  the  effect  of  reducing  computation  time 
by  nearly  a  factor  of  three. 

The  required  bandwidth  increases  with  the  level  of  heterogeneity  in  the  model  (see 
Table  1).  For  model  D  in  figure  3,  a  bandwith  of  ±5  modes  around  the  target  mode  is 


sufficient  for  both  Lg  and  Sn  modes.  However  to  account  for  the  interaction  between 
modes  for  model  E  it  is  desirable  to  have  at  least  a  bandwidth  of  ±7  modes  included  in  the 
calculation.  For  the  highest  level  of  heterogeneity  illustrated  in  fig  3  (model  F)  a  bandwidth 
of  ±9  modes  is  definitely  insufficient  for  full  accuracy.  For  very  heterogeneous  models 
there  is  also  the  possibility  of  high  angle  propagation  effects  and  so  the  size  of  the  mode  set 
has  to  be  increased.  Thus  although  the  coupled  mode  approach  is  not  confined  to  low  level 
heterogeneity,  it  is  at  is  most  effective  when  the  structure  does  not  deviate  too  far  from  the 
reference.  A  convenient  working  limit  would  seem  to  be  about  2  per  cent  deviation  as  in 
model  E. 

In  order  to  understand  the  effect  of  propagation  through  a  heterogeneous  model  on  the 
nature  of  regional  phases  we  need  to  understand  the  characteristics  of  the  surface  wave 
modes  in  the  reference  model.  In  figure  5  we  display  the  group  slowness  behaviour  as  a 
function  of  frequency  up  to  2  Hz  for  the  first  40  Love  modes  on  the  ARANDA  reference 
structure.  The  advantage  of  displaying  the  group  slowness  (the  reciprocal  of  the  group 
velocity)  is  that  the  various  arrivals  can  be  recognised  in  the  time  relationship  that  they 
would  have  on  a  seismic  record. 

The  Sn  phase  can  be  recognised  in  figure  5  as  the  superposition  of  many  slowness 
maxima  and  minima  at  a  slowness  around  0.22  s/km  (4.5  km/s).  All  40  modes  contribute 
to  this  tightly  defined  band  of  slownesses  which  clearly  separates  from  the  other  classes  of 
arrival.  The  onset  of  the  Lg  phase  at  a  slowness  of  0.286  s/km  (3.5  km/s)  is  associated 
with  the  superposition  of  a  number  of  relatively  broad  slowness  minima  arising  from 
relatively  low  order  modes.  Following  this  onset  is  a  more  tangled  skein  of  maxima  and 
minima  with  slownesses  less  than  0.304  s/km  (i.e.  group  velocities  greater  than  3.3  km/s), 
associated  with  about  four  modes  at  each  frequency,  which  will  provide  the  main  Lg 
arrival.  Somewhat  later,  we  have  a  sequence  of  sharply  defined  slowness  maxima 
associated  with  each  mode  in  turn  as  the  frequency  increases.  These  maxima  represent  the 
Airy  phase  for  each  mode  and  will  be  significant  contributors  to  the  Lg  coda:  they  arise 
physically  from  multiple  reflections  at  angles  to  the  vertical  just  greater  than  the  critical 
angle  for  the  crust-mantle  boundary. 

Figure  5:  Group  slowness  beviour  as  a  function  of  frequency  for  the  first  40  Love  modes  of  the  ARANDA 

reference  model  illustrated  in  figure  1. 
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From  the  modal  eigenfunction  patterns  displayed  in  figure  1,  we  can  recognise  that 
those  modes  (3-6)  which  contribute  most  to  the  onset  of  Lg  at  1.5  Hz  have  the  bulk  of  then- 
energy  confined  to  the  middle  and  upper  crust.  The  main  Lg  arrivals  come  from  modes 
7-10  which  sample  the  whole  crust.  The  Airy  phases  for  modes  14,  15  arise  just  before 
the  transition  to  energy  transport  in  the  mantle. 

The  effect  of  varying  the  heterogeneity  level  in  the  structural  models  can  be  well 
illustrated  by  comparing  the  energy  distribution  across  the  modal  sequence  for  a  single 
incident  mode.  In  Table  1  we  show  this  energy  distribution  for  a  set  of  incident  mode  at  0 
km,  at  a  frequency  of  1.5  Hz,  chosen  to  represent  different  parts  of  the  Lg  wave  phase. 
We  consider  horizontal  transmission  through  1000  km  of  structure  for  four  different 
models  C,  D,  E  and  F  with  increasing  levels  of  heterogeneity.  Model  C  has  1  per  cent  of 
hetrogeneity  confined  to  the  crust.  Models  D,  E  and  F  are  illustrated  in  figure  3,  and  have 
the  same  heterogeneity  model  in  the  mantle  with  about  1  per  cent  variation  from  the 
reference  model  (ARANDA,  fig  1),  but  increasing  levels  of  crustal  heterogeneity  (±1  per 
cent  for  D,  ±2  per  cent  for  E,  and  ±5  per  cent  for  F)- 

At  1.5  Hz,  mode  4  is  a  principal  contributor  to  the  onset  of  the  Lg  wave  train.  From 
Table  1  we  see  that  for  the  incidence  of  this  single  mode  at  0  km,  the  majority  of  the  energy 
is  carried  in  the  original  mode  over  the  full"  1000  km  propagation  path  for  the  structures 
with  up  to  2  per  cent  variation  away  from  the  reference  model  (C,D,E).  However,  the 
proportion  of  energy  in  the  original  mode  decreases  as  the  crustal  heterogeneity  increases. 
As  expected  the  mantle  heterogeneity  has  very  little  influence  on  this  group  of  interacting 
modes  which  is  principally  confined  to  the  crust  Once  we  reach  the  highest  level  of 
heterogeneity  (model  F)  the  pattern  of  the  energy  distribution  is  markedly  changed,  the 
energy  maximum  is  shifted  over  two  modes  to  mode  6  which  has  a  slightly  different  group 
velocity  and  energy  is  spread  widely  across  the  whole  suite  of  crustal  modes.  There  is  very 
little  interaction  with  the  fundamental  or  first  higher  modes  because,  as  can  be  seen  from 
figure  1,  their  energy  is  confined  to  the  near  surface  and  so  is  only  sensitive  to  a  small  part 
of  the  total  crustal  heterogeneity. 

For  a  mode  in  the  main  Lg  wave  arrivals  (mode  9),  we  see  that  for  the  two  models  with 
relatively  modest  levels  of  heterogeneity  (C,D)  the  energy  is  concentrated  over  three 
neighbouring  modes  centered  on  the  incident  mode.  These  three  modes  have  group 
velocities  varying  by  about  0.15  km/s  and  so  the  effect  of  the  this  level  of  heterogeneity 
will  be  to  give  a  more  diffuse  maximum  in  the  Lg  wave  amplitude.  The  presence  of  mantle 
heterogeneity  now  has  a  slight  effect  since  it  is  possible  to  couple  into  modes  with  some 
energy  penetration  into  the  uppermost  mantle.  For  the  higher  levels  of  crustal  heterogeneity 
the  incident  energy  gets  spread  out  over  a  broader  range  of  modes  (6  are  significant  for 
model  E  and  at  least  10  for  model  F)  and  the  effect  on  the  Lg  wavetrain  will  be  more 
profound.  For  modest  levels  of  heterogeneity  (around  1  per  cent)  the  best  measure  of  Lg 
wave  strength  will  be  the  integrated  energy  in  a  time  window  spanning  group  velocities  of 
about  3. 6-3. 3  km/s  for  a  typied  continental  situation. 

Similar  results  arise  for  the  Lg  wave  coda  (incident  mode  14)  though  this  case,  uhere 
the  incident  mode  is  most  sensitive  to  lower  crustal  structure,  is  not  as  strongly  affected  by 
the  heterogeneity  as  for  mode  9.  Even  so,  with  increasing  heterogeneity  levels  the 
representation  of  the  propagation  via  the  modes  of  the  reference  structure  requires  extensive 
mode  coupling.  For  the  higher  heterogeneity  levels  there  is  significant  energy  shift  into  Sn 
type  modes,  which  we  can  explain  physically  as  occurring  from  local  changes  in  the  critical 
angle  at  the  crust-mantle  boundary.  Modes  like  14  travel  close  to  the  critical  angle  in  the 
reference  model  and  so  a  slight  change  of  conditions  can  lead  easily  to  energy  transmission 
into  the  mantle. 
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Table  1:  Energy  distribution  across  the  modes  of  the  ARANDA  model  at  1.5  Hz,  after 
transmission  through  1000  km  of  heterogeneous  structure 


a)  Incident  Mode  4  (onset  of  Lg) 
Relative  Mode  Number 
Model:  C 

D 
E 
F 


-3 

-2 

-1 

0 

1 

2 

3 

4 

.000 

.032 

.048 

.689 

.109 

.078 

.003 

.023 

.000 

.032 

.044 

.689 

.109 

.078 

.006 

.023 

m 

.000 

.040 

.212 

.563 

.032 

.036 

.084 

.002 

.  m 

.001 

.098 

.082 

.201 

.057 

.408 

.071 

.024 

b)  Incident  Mode  9  (main  Lg  arrival) 
Relative  Mode  Number:  -4  -3  -2 

Model:  C  .006  .014  .00* 

D  .006  .017  .00i 

E  .008  .023  .14 

F  .078  .102  .00: 


.006  .014  .006  .336  .360  .240  .010  .000  .000 

.006  .017  .006  .336  .348  .250  .012  .000  .000 

.008  .023  .144  .090  .423  .160  .078  .004  .000 

.078  .102  .003  .303  .260  .053  .068  .048  .026 


c)  Incident  Mode  14  (coda  of  Lg) 
Relative  Mode  Number  -4 


-3  -2 


Model:  C 

.000 

.000 

.001 

.032 

.846  .090 

.001 

.000 

.000 

D 

.000 

.000 

.001 

.029 

.828 

.109 

.001 

.000 

.000 

E 

.000 

.000 

.020 

.014 

.397  .476 

.073 

.001 

.005 

F 

.001 

.008 

.006 

.084 

.672 

.109 

.006 

.020 

.014 
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For  incident  modes  4  and  9,  the  nature  of  the  energy  spread  induced  across  the  modes 
of  the  reference  structure  has  been  to  principally  couple  Lg  type  modes  together  and  so  the 
pattern  of  energy  spread  shows  a  skew  towards  higher  mode  numbers  for  mode  4  and 
lower  numbers  for  mode  9.  Whereas,  for  incident  mode  14,  the  coupling  is  largely  into  a 
the  immediate  neighbouring  modes  with  similar  character  and  also  into  the  Sn  males. 

We  should  note  that  even  for  the  moderately  heterogeneous  model  D  (fig  3)  there  has 
been  a  significant  modification  of  the  distribution  of  energy  between  modes  after  passage 
through  1000  km  of  structure.  Such  modifications  over  a  band  of  frequencies  will  lead  to 
theoretical  seismograms  which  will  differ  significantly  from  the  predictions  for  the 
horizontally  stratified  reference  model. 

Studies  of  Three-dimensional  Heterogeneity 

In  addition  to  the  work  on  two-dimensional  structures,  efforts  have  been  made  over  the 
past  year  to  extend  the  mode  coupling  techniques  to  three-dimensional  models.  For 
relatively  weak  heterogeneity  it  has  proved  possible  to  construct  a  approximate 
development  which  should  be  able  to  give  a  good  represenation  of  the  Lg  wave  field  when 
most  significant  energy  travels  close  to  a  direct  path  between  source  and  receiver,  and  so 
multipathing  is  not  too  important.  This  approach  will  require  the  solution  of  partial 
differential  equations  for  the  coupled  modes  in  a  spatial  swath  about  the  direct  path  and 
both  Love  and  Rayleigh  wave  modes  will  need  to  be  considered  at  the  same  time  (cf 
Snieder  1986).  Thus  even  at  1  Hz  at  least  30  coupled  partial  differential  equations  will  be 
required.  The  number  of  modes  needed  at  each  frequency  will  be  at  least  double  that  for 
the  two-dimensional  case  and  also  now  calculations  will  need  to  be  carried  for  a  number  of 
offset  points  in  the  transverse  direction  to  the  path.  This  requirements  will  place  heavy 
demands  on  computational  resources  in  future  work  and  as  a  result  considerable  effort  is 
being  expended  on  the  design  of  suitable  algorithms  for  the  calculations. 


DISCUSSION 

The  current  research  effort  has  been  directed  towards  determining  the  numerical 
characteristics  of  the  Kennett  (1984)  mode-coupling  scheme  for  guided  wave  propagation 
over  extended  distances  at  frequencies  of  1  Hz  and  above.  The  existing  models  are 
two-dimensional  and  calculations  are  carried  out  for  one  frequency  at  a  time. 

The  results  of  the  studies  of  heterogeneous  crustal  models  show  that  reflection  from 
distributed  heterogeneity  is  not  very  important.  However,  both  reflection  and  transmission 
processes  should  be  included,  if  at  all  possible,  in  the  numerical  calculations  in  order  to 
ensure  the  accuracy  of  estimates  of  modal  transmission.  The  coupled  mode  technique 
works  well  with  heterogeneous  models  which  differ  from  a  stratified  reference  model- by  up 
to  2  per  cent  with  a  limited  bandwidth  of  interaction  between  modes.  The  approach  can  be 
used  for  higher  levels  of  heterogeneity  but  a  large  number  of  modes  needs  to  be  considered 
with  consequent  high  computation  costs. 

One  of  the  major  effects  of  crustal  heterogeneity  is  to  introduce  the  possibility  of 
smearing  out  the  main  amplitude  peak  in  the  Lg  wave  train  over  a  band  of  group  velocities. 
As  a  result,  an  effective  measure  of  the  energy  content  of  the  Lg  waves  will  be  to  consider 
the  integrated  amplitude  along  the  traces  beween  group  velocities  of  3.6  and  3.3  km/s.  The 
influence  of  mantle  heterogeneity  on  Lg  wave  propagation  was  found  to  be  quite  small. 
However,  any  conversion  into  Sn  type  modes  will  give  arrivals  which  will  have  apparent 
group  velocities  of  between  4.2  and  3.7  km/s,  depending  on  the  position  where  conversion 
occurs.  Such  arrivals  will  appear  in  a  portion  of  the  seismic  record  which  would  be 
predicted  to  be  very  quiet  in  stratified  medium  calculations,  and  as  a  result  may  appear  to  be 
more  prominent  than  expected. 


The  costs  of  performing  the  coupled  mode  calculations,  at  fixed  step  length  in  the 
horizontal  direction,  would  rise  as  the  square  of  the  number  of  modes  employed  (and 
therefore  approximately  as  the  frequency)  if  full  modal  coupling  is  allowed.  The 
introduction  of  a  limited  bandwidth  of  interaction  reduces  the  cost  but  still  it  will  be  difficult 
to  push  the  calculations  to  very  high  frequency.  At  2  Hz,  about  35  modes  are  needed  for  a 
typical  continental  model  to  give  an  adequate  representation  of  the  effects  of  heterogeneity 
on  either  Love  or  Rayleigh  waves,  and  this  would  increase  to  about  70  modes  for  4  Hz 
propagation.  As  the  wavelength  reduces  the  horizontal  step  length  in  the  solution  of  the 
differential  equations  needs  to  be  reduced  to  maintain  accuracy.  This  gives  a  further 
increase  in  computational  cost  approximately  linear  in  frequency.  The  numerical  algorithms 
used  for  the  solution  of  the  matrix  Ricatti  equations  for  reflection  and  transmission  work 
well  with  large  numbers  of  modes,  but  the  differential  equations  are  non-linear  and  there  is 
the  possibility  of  loss  of  precision  if  very  large  numbers  of  coupled  equations  are 
considered. 

Nevertheless  the  coupled  mode  approach  to  the  calculation  of  guided  wave  propagation 
has  the  potential  to  enable  the  character  of  Lg  wave  propagation  to  be  determined  for 
classes  of  model  which  cannot  easily  be  treated  by  any  other  means.  Both  deterministic  and 
statistical  models  of  heterogeneity  can  readily  be  incorporated  into  the  computation  scheme. 
The  coupled  mode  scheme  does  therefore  allow  an  assessment  of  structural  effects  on 
regional  phase  propagation  which  is  complementary  to  existing  techniques  using 
horizontally  stratified  models. 

The  next  phase  of  work  will  be  directed  towards  using  the  experience  from  the 
two-dimensional  models  to  design  a  computational  scheme  for  mode-coupling  in 
three-dimensions  for  both  Love  and  Rayleigh  waves.  This  scheme  will  take  advantage  of 
approximations  such  as  limited  bandwidth  of  interaction  to  reduce  total  computation  time. 
We  will  also  pursue  the  problem  of  generating  effective  and  consistent  theoretical 
seismograms  by  combining  mode-coupling  results  from  a  number  of  frequencies. 
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APPENDIX:  Modal  evolution  aquations 
1 .  Two— dimension*!  problmms 

We  present  here  a  oriet  derivation  of  the  coupled  mode 
technique  introduced  by  Kennett  (1984). 

We  will  •firstly  restrict  attention  to  a  two-dimensional  iy 
heterogeneous  structure  and  work  with  a  fully  anisotropic  medium 
for  which  a  very  compact  notation  is  available.  The  heterogeneity 
is  assumed  to  be  invariant  in  the  transverse  direction  to  the  wave 
vector  A  ,  and  so  the  direction  of  this  vector  will  remain 

constant  even  in  the  anisotropic  case  where  the  ray  direction 
associated  with  a  particular  mode  does  not  lie  along  A  . 

We  consider  the  heterogeneity  as  a  deviation  from  the 
properties  of  a  reference  model.  This  perturbation  is  not 
necessarily  small  but  must  be  bounded  in  space.  We  seek  a 
representation  of  the  displacement  field  in  terms  of  a 
superposition  of  the  modal  eigenfunctions  of  the  reference  model 
with  modal  coefficients  which  vary  with  position. 

We  work  in  a  cartesian  coordinate  system,  with  propagation 
in  the  Jt  -  z,  plane  and  introduce  the  displacement  vector 

%  =  (  k/j, ,  k.y,b/z)  .  We  also  need  to  specify  the  traction  field. 

In  the  case  of  a  modal  wavetrain  propagating  in  the  JC  —direction 
we  are  concerned  with  continuity  of  traction  on  vertical  planes 
and  so  we  take  £t  •  ( T„  ,  r,z.  t„)  -  i  n  addition  the  material 

continuity  requirements  mean  that  we  also  need  to  be  able  to  refer 
to  the  traction  on  a  horizontal  plane  dj  =•  ('fj,  ,  TJt .  Tu )  -  For 

the  2— D  situation  we  need  .?  number  of  combinations  of  elastic 
moduli  which  may  be  represented  compactly  in  terms  of  matrices 
C«v  such  that 


(Cy)ll  =  C*UJ 


where  the  are  the  anisotropic  moduli. 

The  equations  of  motion  and  the  stress-strain  conditions  can 
be  cast  into  a  form  where  the  derivatives  with  respect  to  the 
propagation  direction  JC  appear  only  on  the  left  hand  side  of 
the  equations 

a  l  tr\  /  A~t  \  !  z\  ti\ 


in  the  presence  of  a  volume  force  contribution  / 

The  differential  operators  A..  do  not  depend  on  the  horizontal 
derivatives  of  the  material  properties  and  for  the  2-D  case  nave 
the  explicit  form 

/  \  /  -  C„  C,SA&  C„  ) 

(  A#  J  (  1  -  4  £$33^  C"  ' 

where  we  have  written  for  A/?z.  and  set 

Qs  *  CSt  '  C  Cn  C,i  (3) 

The  unclosed  brackets  in  (2)  indicate  that  the  operator  acts  to 
its  right.  Across  an  interface  Z  -  ,  the  Vt)  component  of 


13 


will  be  continuous,  but  7n  ,  7rx 
continuous.  On  the  other  hand 
be  represented  in  terms  of  if, 


are  not  required  to  be 

will  be  continuous  and  can 
and  as 


Q..  4  z  **  cJt  c;‘£ 


The  equations  <l)-<3)  apply  in  an  arbitrarily  heterogeneous  medium 
but  by  themselves  provide  no  direct  link  to  the  surface  wave  case 
we  wish  to  consider. 

When  the  heterogeneous  medium  does  not  deviate  too  strongly 
from  the  reference  medium,  we  can  envisage  that  it  may  be  possible 
to  find  a  representation  in  terms  of  the  surface  wave  modes  of  the 
reference  structure.  We  here  consider  the  case  of  a  fixed, 
stratified,  reference  medium  as  used  in  the  examples  of  Kennett 
(1984)  and  Kennett  &  Mykkeltveit  (1984). 

At  each  position  Jc  we  consider  cutting  the  actual  structure 
along  a  vertical  plane  and  then  weld  on  the  reference  structure  to 
ensure  continuity  of  displacement  jf  and  traction  .  In  this 

reference  structure  we  can  now  represent  the  displacement  and 
traction  at  x  as  a  superposition  of  modal  contributions  and 
write 


(  )  *  Z  CM  e*/>(i£rx)  (* 


(  K.*) 


where  the  ff°  are  the  displacement  eigenfunctions  for  the 
reference  medium  and  the  horizontal  tractions  are  derived 

from  *K-“  .  The  sum  is  to  be  taken  over  all  relevant  modes  at 

frequency  co  .  In  order  to  give  a  full  representation  of  the 
displacement  this  sum  will  normally  involve  both  forward  and 
backward  travelling  waves  (i.e.  both  positive  and  negative  ). 

As  discussed  by  Kennett  (1984)  we  can  arrange  the  deepest  part 
of  the  model  of  the  reference  medium  so  that  (5)  is  a  sum  over 
an  infinite  set  of  orthogonal  modes,  and  achieve  a  complete 
representation  of  the  seismic  wavefield.  For  example,  the 
constructive  interference  of  sufficient  numbers  of  surface  wave 
modes  will  synthesise  body  wave  phases  (see  e.g.  the  examples  in 
chapter  1 1  of  Kennett  (1983)).  However,  especially  at  higher 
frequencies,  this  requires  a  very  large  number  of  modes.  In  the 
varying  medium  we  wish  to  concentrate  on  the  surface  wave  field 
and  so  would  like  to  be  able  to  work  in  terms  of  a  restricted  (and 
not  too  large)  set  of  modes.  Such  a  truncation  of  the  expansion 
(S)  imposes  limitations  on  the  size  and  character  of  the 
heterogeneity  which  can  be  tackled. 

We  recall  that  T„  ,  TtM  are,  in  general,  not  continuous' across 
horizontal  interfaces  and  so  if  the  material  discontinuities  in 
the  heterogeneous  model  do  not  coincide  with  those  in  the 
reference  model,  we  have  the  problem  of  representing  a  jump  in 
■fr  across  an  interface  by  a  superposition  of  continuous 
traction  vectors  .  A  satisfactory  fit  <-an  be  achieved  with 

many  modes,  but  the  accuracy  of  the  representation  may  be  limited 
with  a  restricted  mode  set  (cf.  Fourier  series).  As  a  result  it 
is  desirable  that  the  major  discontinuities  in  both  the 
heterogeneous  and  reference  models  should  be  coincident. 
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In  the  reference  structure  the  modal  contr 1 buti ons  (5) 
satisfy  the  coupled  equations 


zm  { 

y  \  t‘ 

\  Zrr  / 


-  2  4-  {•'  -  <>, 


(6) 


in  the  absence  of  any  applied  force,  with  coefficients 
Cr  £je/>  .  We  have  written  A° 


Ct- 


to  indicate  the  form  of  the  differential  operators  for  the 
reference  medium. 

In  a  laterally  varying  medium  we  have  to  require  the  modal 
coefficients  to  vary  with  position.  On  expanding  the 
differential  operators  A  as  A 0  *■  A  A  ,  to  separate  off  the 
reference  medium  contribution,  we  find 


r/  I  a ;/  r  [tj 


<  7 ) 


*  m  Z  4 


(8) 


in  which  we  can  recognise  the  left  hand  side  as  having  the  same 
form  as  in  (6),  but  we  have  in  addition,  on  the  right  side,  a  term 
which  is  equivalent  to  a  generalised  volume  force  applied  to  the 
reference  structure. 

As  well  as  the  equation  of  motion  we  require  the  wavefield  to 
satisfy  certain  continuity  conditions  at  internal  and  external 
boundaries  where  the  elastic  moduli  are  discontinuous.  The  most 
general  case  of  an  interface  in  the  heterogeneous  medium  will  be 
titled  and  displaced  from  those  of  the  reference  medium.  Continuity 
of  displacement  is  assured  at  any  interface  from  the  continuity 
properties  of  the  displacement  eigenfunctions.  We  also  require 
that  the  tractions  normal  to  any  interface  should  be  continuous. 

We  will  assume  that  all  such  surfaces  do  not  vary  rapidly  in  the 
horizontal  direction.  Then  tor  an  interface  in  the 

laterally  varying  medium  described  by  the  function  A.{x)  ,  we  may 

work  to  first  order  in  the  slope  (  A,  )  ,  and  require 


[tii  *  cItc;;)Z  c,u>i;r  .  q„z  zrw  z  »■,•]•  ° 


(9) 


where  I  denotes  the  identity  matrix  and  the  square  brackets 
indicate  the  jump  in  the  enclosed  quantities  across  the  interface, 
evaluated  from  bottom  to  top. 

If  we  now  separate  out  the  contribution  from  the  properties 
of  the  reference  medium  from  the  local  variations,  we  may  write 


[(C,C-fZcrMt-  y  <2 ;zty*>is;l 

r  rr  r 


-//■ 


(  lO) 


The  difference  terms  such  as  d  (Cj,  C,  )  in  (10)  indicate  the 
discrepancy  between  the  combinations  of  elastic  moduli  for  the 
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heterogeneous  and  re+erence  models.  In  the  absence  of  any  lateral 
heterogenei ty  the  expression  on  the  left  hand  side  of  (10)  would 
vanish.  The  presence  of  the  heterogeneity  is  thus  equivalent  to 
inducing  a  traction  discontinuity  into  the  reference  medium,  along 
the  line  of  the  interface.  Such  a  traction  discontinuity  is 
equivalent  to  a  localised  volume  force  along  the  interface 


fr]  S(z-a,m) 


(see  e.g.  section  3.1  of  Aki  &  Richards  (1980)). 

We  can  therefore  represent  the  combination  of  the  equations  of 
motion  and  the  boundary  conditions  in  the  laterally  varying 
medium  by  means  of  the  fallowing  system  of  equations 


Z  4 


* 


.  (11) 


where  the  force  term  is  a  summation  of  interface  contributions 

/„  ■  £  [  u  (q  a ,  4  //  /;  *  ^<3, 4  !■;  I  *» 

With  the  abstraction  of  this  force  contribution  the  interface 
conditions  on  the  wavefield  are  reduced  to 

l  Z  O*  K‘  l  ■  o.  [ Z  X>f  (4  cp-£  ,  Qigrll  - 

where  the  subscript  -n  indicates  the  jump  at  any  interface. 

These  constraints  will  be  automatically  satisfied  by  the 
representation  (5).  Equation  (11)  can  be  further  simplified 
because  the  modal  eigenfunctions  are  just  the  free  solutions  for 
the  reference  structure  and  so  the  left  hand  side  of  the  equation 


vanishes. 

This  leaves 
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We  now  exploit  the  orthogonality  relation  between  different 
modal  eigenfunctions  for  the  reference  medium 

*  /  f  ~  iy,-  (13) 

to  get  a  set  of  coupled  first  order  differential  equations  for  the 
modal  coefficients  f  <Z-  f  ,  appropriate  for  both  elastic  and 
inelastic  heterogeneity.  These  equations  can  be  written  in  the 
form 

£  •■£  <?,  Ja  f  gj.  -  4(c~cM<l 

-  /-  ‘  JCZ  C>£W 

“Z  4-  UX  -  '“Li* '  4(4001  ■  <14’ 

where  we  have  written 

(~  kf.  *■)  ■ 
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Equation  (14)  can  be  simplified  by  integration  by  parts,  to  yield 

-iZ  £r  /'* 

-.4’  ^  Z-'Oix;  -  *  A&OKll 
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4  -f  4  f"/.  CL 


(15) 


In  which  we  see  that  there  is  no  explicit  interface  term  in  the 
case  of.  a  flat  boundary  (  •  0  ).  However,  the  modification 

of  an  interface  will  appear  through  the  integral  term  in  (15), 
where  the  properties  of  the  laterally  varying  medium  differ  from 
the  reference  and  there  will  be  a  significant  contribu-ion  to  the 
integral  which  will  vary  with  horizontal  position.  In  a 
first-order  perturbation  theory  the  whole  of  this  effect  would  be 
projected  onto  a  specific  interface  term,  but  here  we  are  able  to 
make  a  more  detailed  allowance  for  the  behaviour. 

The  combination  of  the  modal  representation  (S)  with  the 
imposition  of  continuity  of  traction  at  each  interface  thus 
leads  to  the  set  of  coupled  equations 


<?x 


1  ■  ‘  2  V  cx/  /-  x)  ax/  f-itr*)  4- 


(16) 


for  the  modal  coefficients,  where  the  coupling  coefficients 
between  the  modes  can  be  written  as 

-4,  -  /"^  / 1,  j*: &  * 

- z -  <,'■<>  sw  4*  f 
'  ?  4.  /  Ft-  4-X 


(17) 

The  coupled  first-order  equations  (16)  are  not  very  easy  to 
solve  because  we  have  a  two— point  boundary  value  problem  with  both 
reflected  and  transmitted  modes  to  be  determined.  Kennett  (1984, 
section  3)  has  shown  that  an  effective  procedure  is  to  work  with 
the  reflection  and  transmission  matrices  tor  a  sequence  of  models 
encompassing  increasing  portions  of  the  heterogeneous  medium. 

This  leads  to  an  initial  value  problem  for  two  coupled  matrix 
Ricatti  equations.  These  nonlinear  differential  equations  for  the 
reflection  and  transmission  matrices  can  be  readily  solved 
numerical 1 y. 


2  T hr »m-d immn s ion ml  Problems 

In  a  medium  whose  properties  vary  in  three  dimensions, 
it  is  possible  to  generate  evolution  equations  -for  the 
displacement  and  horizontal  tractions  with  respect  to  the 
horizontal  coordinates.  The  analysis  parallels  the  two-dimensional 
case,  but  now  one  has  to  take  account  of  variation,  in  both  the 
material  properties  and  the  wave-field,  transverse  to  the  sirection 
o-f  current  interest. 

If  we  concentrate  on  the  evolution  m  the  -X  -direction  we  -find 


Cm" 

i  -  j>z  [ 

-SzWuSy  -  <?y  ~  >4 /£, C 


where 


4  -  4  C  - 


4,  • 


In  equation  (18),  the  terms  have  been  written  in  an  order  which 
emphasises  those  contributions  which  appear  in  the  two-dimensional 
-formulation.  Previously  we  were  able  to  eliminate  the  dependence  on 
depth  z  and  convert  from  evolution  equations  in  displacement  and 
traction  to  equations  in  the  modal  expansion  coefficients  by 
making  use  of  the  modal  orthogonal i ty  properties  for  modes 
travelling,  in  the  same  (or  opposite)  directions.  With  3-D 
heterogeneity  we  now  have  the  possibility  of  modes  travelling  in 
different  directions  and  also  coupling  between  Love  and  Rayleigh 
modes.  Unf ortunately,  we  cannot  write  down  any  simple  form  for 
the  representation  of  the  wavefield  as  a  sum  of  modal 
contributions. 

If  we  consider  a  situation  where  the  heterogeneity  is  not  too 
strong  and  concentrate  attention  on  propagation  dominantly  along 
the  jc  direction,  we  can  generate  a  useful  approximation  by 
neglecting  the  angular  dependence  of  modal  interaction.  By  this 
means,  we  obtain  a  generalisation  of  the  2-D  equations  with  a  weak 
dependence  on  the  transverse  (  jr  )  direction.  The  partial 
differential  equations  in  x  for  the  .modal  coefficients  are  then  to 
be  integrated  over  a  swath  surrounding  the  direct  propagation  path 
from  source  to  receiver. 


[VftlAi 


» 


DISTRIBUTION  OF  RESEARCH  RESULTS 


I 


a)  PUBLICATIONS: 

Lg-wave  propagation  in  a  heterogeneous  crust, 

(B.L.N.  Kennett)  In  preparation 

for  submission  to  the  Bulletin  of  the  Seismological  Society  of  America 


b)  PRESENTATIONS: 

On  representing  surface  wave  propagation  in  laterally  varying  media  by  truncated  modal 
expansions 

DARPA/AFGL  Seismic  Research  Symposium 
Nantucket,  Massachusetts,  June  1987 

Guided  waves  in  laterally  heterogeneous  media 

Seminar,  Hawaii  Institute  of  Geophysics,  August  1987 


•iTiMVOiII 


»Vi\ 


CONTRACTORS  (United  Scat. 


Professor  Keiiti  Aki 

Center  Cor  Earth  Sciences 

University  of  Southern  California 

University  Park 

Los  Angeles,  CA  90089-0741 

Professor  Charles  B.  Archambeau 
Cooperative  Institute  Cor  Resch 
in  Environmental  Sciences 
University  of  Colorado 
Boulder,  CO  80309 

Dr.  Thomas  C.  Bache  Jr. 

Science  Applications  Int'l  Corp. 

102L0  Campus  Point  Drive 

San  Diego,  CA  92121  (2  copies) 

Dr.  Douglas  R.  Baumgardt 
Signal  Analysis  4  Systems  Div. 
ENSC0,  Inc. 

5400  Port  Royal  Road 
Springfield,  VA  22151-2388 

Dr.  S.  Bratt 

Science  Applications  Int'l  Corp. 
10210  Campus  Point  Drive 
San  Diego,  CA  92121 

Dr.  Lawrence  J.  Burdick 
Woodward-Clyde  Consultants 
P.0.  Box  93245 

Pasadena,  CA  91109-3245  (2  copies) 

Professor  Robert  W.  Clayton 
Seismological  Laboratory /Div .  of 
Geological  5.  Planetary  Sciences 
California  Institute  of  Technology 
Pasadena,  CA  91125 

Dr.  Vernon  F.  Cormier 

Department  of  Geology  4  Ceophysics 

U-45 ,  Roun  207 

The  University  of  Conneticut 
Storrs,  Connecticut  06268 

Dr.  Zoltan  A.  Der 
ENSC0,  Inc. 

5400  Port  Royal  Road 
Springfield,  VA  22151-2388 

Professor  John  Ferguson 
Center  for  Lithospheric  Studies 
The  University  of  Texas  at  Dallas 
P.0.  Box  830688 
Richardson,  TX  75083-0688 


v.iuiv  XT 

Professor  Stanley  Flatce' 

Applied  Sciences  building 
University  of  California,  Santa  Cruz 
Santa  Cruz,  CA  95064 

Professor  Steven  Crand 
Department  of  Ceology 
245  Natural  History  Building 
1301  West  Green  Street 
Urbana,  IL  61801 

Professor  Roy  Greenfield 
Geosciences  Department 
403  Deike  building 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Professor  David  G.  Harkrider 
Setsmological  Laboratory 
Div  of  Geological  a  Planetary  Sciences 
California  Institute  of  Technology 
Pasadena,  CA  91125 

Professor  Donald  V.  Helmberger 
Seismological  Laboratory  ^ 

Div  of  Geological  6  Planetary  Sciences 
California  institute  of  Technology 
Pasadena,  CA  91125 

Professor  Eugene  Herrin 
Institute  for  the  Study  of  Earth 
a  Man/Ceophysical  Laboratory 
Southern  Methodist  University 
Dallas,  TX  75275 

Professor  Robert  B.  Herrmann 
Department  of  Earth  &  Atmospheric 
Sciences 

Saint  Louis  University 
Saint  Louis,  M0  63156 

Professor  Lane  R.  Johnson 
Setsmographic  Station 
University  of  California 
Berkeley,  CA  94720 

Professor  Thomas  H.  Jordan 
Department  of  Earth,  Atmospheric 
and  Planetary  Sciences 
Mass  Institute  of  Technology 
Cambridge,  MA  02139 

Dr.  Alan  Kafka 
Department  of  Ceology  & 

Ceophys  tes 
Boston  College 
Chestnut  Hill,  MA  02167 


'iVWiJ 


'V.n'.nVA’ViVaIWai 


^  Min  m 


m,  w  v,vv.>;  >/  ^  >V.  A  BC  TOVT ^  W  7  TJ  W  WM 


Professor  Leoa  Knoporf 
University  of  California 
Inscicute  of  Geophysics 
£.  Planetary  Physics 
Los  Angeles,  CA  90024 

Professor  Charles  A.  Langston 
Geosciences  Department 
403  Deike  Building 
The  Pennsylvania  Scate  University 
University  Park,  PA  16802 

Professor  Thorne  Lay 
Department  or  Geological  Sciences 
LuOo  C-G.  Little  Building 
University  of  Micnigun 
Ann  Harbor,  Hi  w8109-106J 

Dr.  RnndoLph  Martin  III 
Mow  Lug  land  Research,  Inc . 

P.0.  Box  857 
Norwich,  VT  05055 

Dr.  Gary  McCarter 
Mission  Research  Corp. 

735  Scate  Street 
P . 0 .  Drawer  719 

Santa  Barbara,  CA  93102  (2  copies) 

Professor  Thomas  V.  McEvilly 
Se  i  smog  rapiiic  Station 
University  of  California 
Berkeley,  CA  94720 

Dr.  Keith  L.  McLaughlin 
S-GUBGD, 

A  Division  of  Maxwell  Laboratory 
P.0.  Box  1620 
La  Jolla,  CA  92038-U20 

Professor  William  Menke 
Lumont-Doherty  Ceological  Observatory 
of  Columbia  University 
Palisades,  NY  I09b4 

Professor  Brian  J.  Mitchell 
Department  of  Earth  6  Atmospheric 

Sc  ieilCeS 

Saint  Louis  University 
Saint  Louis,  MO  o3156 


Mr.  Jack  Murphy 
S-CUUED 

A  Division  of  Maxwell  Laboratory 
11800  Sunrise  Valley  Drive 
Suite  L2L2 

Res  ton ,  VA  2209  1  (  2  copies) 


-3- 


v  wv  wvwvwvinraCT 


v»jrry'jiG.-«Q1-»v»yrv 


Professor  J.  A-  Oroutt 

institute  of  Oi_-oi.hysi.es  and  Plane  Cary 
Physics,  .4-205 

Scripps  Institute  of  Oceanography 
Univ.  of  California,  San  Diego 
La  Jolla,  CA  92093 

Professor  Keith  Priestley 
University  of  Nevada 
Mackay  School  of  Mines 
Reno ,  NV  S9  5  57 


Wilmer  Rivers 
Teicdyue  Ceotech 
314  Moncgoiaery  StreeC 
Alexandria,  VA  22314 

Professor  Charles  G.  Sauuuis 
Center  tor  Starch  Sciences 
University  of  Sauthern  California 
University  Park 
uJS  Ange  leS  ,  u  A  90089—0741 

Dr.  Jeffrey  L-  Stevens 
S-CUULD, 

A  Division  of  Maxwell  Laboratory 
P.0.  Box  1620 
La  Jolla,  CA  92038-1620 

Professor  Brian  Stomp 

institute  for  the  Study  of  Earth  4  Man 
CcophysloaL  Laboratory 
Southern  MecnodLst  University 
Julius,  TM  75275 

Professor  Ta-liang  Teng 
CeuCer  for  Earth  Sciences 
University  of  Souchern  California  - 
University  PurK 
Los  Angeles,  Lii  90089—0741 

Professor  M.  Nafi  Toksoz 
Earth  Resources  Lab 
Dept  of  Earth,  Atmospheric  and 
Planetary  Sciences 

Massachusetts  Institute  ot  Technology 
-•2  Curie  ton  Street 
Cambridge,  MA  02142 

Professor  Terry  C.  Wallace 
Department  of  Geosciences 
Building  ;/ 11 
University  of  Arizona 
Tucson,  .4Z  85721 

WeLdlinger  Associates 
ATTN:  Ur.  Gregory  Wojclk 

t/20  Hansen  Way,  Suite  1UU 
Palo  Alco,  CA  941)04 
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Rockwell  Internal: 1 1  Science  Center 
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Professor  Shelton  S.  Alexander 
Consciences  Department 
403  Deike  Building 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Dr.  Ralph  Archuleta 
Department  of  Ceological 
Sciences 

Univ-  of  California  at 
Santa  Barbara 
Sanca  Barbara,  CA 

Dr.  Muawia  Barazangi 
Geological  Sciences 
Cornell  University 
Ithaca,  NY  14853 

J.  Barker 

Department  of  Geological  Sciences 
Stace  University  of  New  York 
at  Binghamton 
Vestal,  NY  13901 

Hr.  William  J.  Best 
907  Westwood  Drive 
Vienna,  VA  22180 

Dr.  N.  Biswas 
Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 

:.>r  v  i .  Bollinger 

» rtuient  of  Geological  Sciences 
■girila  Polytechnical  Institute 
.  .04-.  erring  Hall 
Blackaourg,  VA  24061 
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Dr.  James  Bulau 

Rockwell  Int 1 1  Science  Center 

1049  Carnino  Dos  Rios 

P.0.  Box  1085 

Thousand  Oaks,  CA  91360 

Mr.  Roy  Burger 
1221  Serry  Rd . 

Schenectady,  NY  12309 
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Dt .  Robert  iiur  ridge 
Schlumbe  rcer-Do  i  i  Resell  Ctr. 
uld  Quarry  Road 
Ridgefield,  CT  06877 

Science  Horizons,  Inc. 

ATTN:  Dr.  Theodore  Cherry 
7  L 0  Lncinitas  illvd.,  SuiCe  101 
Cneinicas,  CA  92024  (2  copies) 

Professor  Jon  F.  Claerbout 
Professor  Amos  Nur 
Dope,  of  Ceophysics 
Stanford  University 
Stanford,  CA  94305  (2  copies) 

Dr.  Anton  W.  Dainty 
AL’GL/LUU 

l  la  ns  com  AFli,  AA  0173  L 
Dr.  S  t  e veil  Day 

Dept,  of  CeologicaL  Sciences 
San  Diego  State  U. 

San  Diego,  CA  92182 

Professor  Adam  Dciewonski 
llofiuian  Laboratory 
Harvard  University 
20  Oxford  St. 

Cambridge,  AA  02138 

Professor  John  Ebel 

Dept  of  Geology  a  Geophysics 

boston  College 

dies tuuc  Hill,  MA  02167 

Dr.  Alexander  Florence 
Silt  International 
333  Ra venswood  Avenue 
Menlo  Park,  CA  94025-3493 

Dr.  Donald  Forsyth 
Dept,  of  Ceologieal  Sciences 
arown  University 
Providence,  HI  02912 

Dr.  Anthony  Cangi 
Texas  AaM  University 
Department  of  Geophysics 
Coliege  Station,  TX  77843 

Dr.  Freeman  Gilbert 
Institute  of  Ceophysics  & 
Planetary  Physics 
Univ.  of  California,  San  Diego 
P • 0 .  box  LU9 
La  Jolla,  CA  92037 
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i-ir  .  Edward  Ciller 

Pacific  Seirra  Research  Corp  . 

1401  Wilson  Boulevard 
Arlington,  VA  22209 

Dr .  Jeffrey  W  .  Given 
Sierra  Geophysics 
11255  Kirkland  Way 
Kirkland,  WA  98033 

Dr  .  Henry  L .  Gray 

Associate  Dean  of  Dedman  College 
Department  of  Statistical  Sciences 
Southern  Methodist  University 
Dallas,  TX  75275 

Rong  Song  Jih 
'i'eledyne  Geotech 
3  1-1  Montgomery  Street 
Alexandria ,  Virginia  22314 

Professor  F.K.  Lamb 
University  of  Illinois  at 

Urbana-Champaign 
Department  of  Physics 
1110  West  Green  Street 
Urbuna,  IL  61801 

Dr  .  Arthur  Lerner-Lam 

Lamont-Donerty  Geological  Observatory 
of  Coluitbia  University 
Palisades,  NY  10964 

Dr.  L.  Timothy  Long 
School  of  Geophysical  Sciences 
Georgia  Institute  of  Technology 
Atlanta,  GA  30332 

Dr  .  Peter  Malin 

University  of  California  at  Santa  Barbara 
Institute  for  Central  Studies 
Santa  Barbara,  CA  93106 

Dr.  George  R.  Mellman 
Sierra  Geophysics 
11255  Kirkland  Way 
Kirkland,  WA  98033 

or  .  Bernard  Minster 

Institute  of  Geophysics  and  Planetary 
Physics,  <\-2  05 

Scnpps  Institute  of  Oceanography 
Univ  .  of  California,  San  Diego 
La  Jolla,  CA  92093 

Professor  John  Nabelek 
College  of  Oceanography 
Oregon  State  University 
Corvallis.  OR  97331 
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Dr  .  Gaza  Nagy 
U.  California,  San  Diego 
Dept  of  Aines,  M  .S  .  3-0  10 
La  Jolla,  CA  32093 

Dr  .  Jack  Oliver 
Department  of  Geology 
Cornell  University 
Ithaca,  NY  14850 

Dr.  Robert  Phinney/Dr  .  F.A.  Dahlen 
Dept  of  Geological 

Geophysical  Sci  .  University 
Princeton  University 
Princeton,  NJ  08540  (2  copies) 

RADIX  Systems,  Inc. 

Attn:  Dr.  Jay  Pulli 

2  Taft  Court,  Suite  203 
Rockville,  Maryland  20850 

Professor  Paul  G.  Richards 
Laiuont-Doherty  Geological 
Observatory  of  Columbia  Univ  . 

Palisades,  NY  10964 

Dr  .  Norton  Rimer 
S-CUUED 

A  Division  of  Maxwell  Laboratory 
P  .0  .  1620 

La  Jolla,  CA  92038-1620 

Professor  Larry  J .  Ruff 
Department  of  Geological  Sciences 
1006  C.C.  Little  Building 
University  of  Michigan 
Ann  Arbor,  MI  48109-1063 

Dr.  Alan  S.  Ryall,  Jr. 

Center  of  Seismic  Studies 
1300  North  17th  Street 
Suite  1450 

Arlington,  VA  22209-2308  (4  copies) 

Dr  .  Richard  Sailor 
TASC  Inc . 

55  Walkers  Brook  Drive 
Reading,  MA  01867 

Thomas  J.  Sereno,  Jr. 

Service  Application  Int'l  Core. 

10210  Campus  Point  Drive 
San  Diego,  CA  92121 

Dr  .  David  G  .  Simpson 
Lamont-Doherty  Geological  Observ  . 

of  Columbia  University 
palisades,  NY  10964 
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Dr .  Bob  Smith 
Department  of  Geophysics 
University  of  Utah 
1400  East  2nd  South 
Salt  Lake  City,  UT  84112 

Dr  .  S  .  W  .  Smith 
Geophysics  Program 
University  of  Washington 
Seattle,  WA  98105 

Dr  .  Stewart  Smith 
IRIS  Inc  . 

161o  N.  Fort  Myer  Drive 
Suite  1440 

Arlington,  VA  22209- 

Rondout  Associates 
ATTN:  Dr.  George  Sutton, 

Dr.  Jerry  Carter,  Dr.  Paul  Pomeroy 
P  .0  .  Box  224 

Stone  Ridge,  NY  12484  (4  copies) 

Dr  .  L  .  Sykes 

Lamont  Doherty  Geological  Observ  .  ^ 

Columbia  University 
Palisades,  NY  10964 

Dr  .  Pradeep  Talwani 
Department  of  Geological  Sciences 
University  of  South  Carolina 
Columbia,  SC  29208 

Dr.  R.  0 .  Tittmann 

Rockwell  International  Science  Center 
1049  Camino  Dos  Rios 
P  .0  .  Box  1085 
Thousand  Oaks,  Ci ft,  91360 

Professor  John  H.  Woodhouse 
Hoffman  Laboratory 
Harvard  University 
20  Oxford  St . 

Cambridge,  MA  02138 

Dr  .  Gregory  B .  Young 

ENSCO ,  Inc . 

5400  Port  Royal  Road 
Springfield,  VA  22151-2388 
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Ter conophys cc|ue  -  I . R . I . G . M-b - P • 
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CedeX  1  RANCH 

ur-  ililiuar  bunguni/NTNl'/NORSAR 
P .  0  -  box  j  I 

Norwegian  Council  or  Science, 

Luduscry  and  Research,  NORSAR 
M - ‘d U J 7  Rjeller,  NORWAY 

Dr.  Michel  Cuuipillo 
r.R. L .C.M.-U.P.  60 
jO-'.ud  Sc.  Harcin  D’ Hems 
CedeX,  1  RANCH 

Dr.  kin- Yip  Chun 
Geophysics  DivLcion 
PhySCCS  DepaCCUienC 
UuiveL’iiCy  or  iYi'uiiLi) 
undario,  C ANaDA  M5S  1A7 

Dr  .  .vlan  Douglas 
Miniacry  or  Detense 
.i  iacRues  C  ,  brinupCGu, 

Read  r  nr  RC7— -»KS 
UNITED  .<  INC  DOM 

Dr.  Man! red  Heuger 

i'ed  .  Lnsc.  For  Ceosciences  d  Nac‘1  Res. 
I'u.. cr’ach  5101  S3 
D-3uOO  Hannover  51 
i  CORRAL  RHPUbLIC  OH  CHRMANY 

Dr.  H  .  Iluse  bye 
NTNH/  NORSAR 
P .  0 .  liox  51 

N-d0o7  Rjeller,  NORUAY 
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Ms  .  Eva  Johunnisson 
Senior  Research  Officer 
National  Defense  Research  Inst  . 

P  .0  .  Box  27322 
S-102  5-1  Stockholm 
SWEDEN 

Tornioti  Kvaerna 
N'i'NF/  NORSAR 
P  .0  .  Box  5  1 

N-2007  Kjeller,  NORWAY 

Mr.  Peter  Marshall,  Procurement 
Executive,  Ministry  of  Defense 
Blackriest,  Drxmpton, 

Roadiny  FG7--1RS 
UNITED  KINGDOM  (3  copies) 

Dr  .  Ben  Menaheim 

Ueizman  Institute  of  Science 

Rehovot,  ISRAEL  95  1729 

Dr  .  Svein  Mykkeltveit 
NTNF/HORSAR 
P  .0  .  Box  5 1 

N-2007  Kjeller,  NORWAY  (3  copies) 

Dr  .  Robert  North 
Geophysics  Division 
Ceoloyical  Survey  of  Canada 
1  Observatory  crescent 
Ottawa,  Ontario 
CAWnDA ,  K1 A  0Y3 

Dr  .  Frode  Rinydal 
NTNF/UORSAR 
P  .0  .  Box  5  1 

N-2007  Kjeller,  NORWAY 
Dr  .  Jory  Schlitteuhnrdt 

federal  Inst,  for  Geosciences  a  Nat’l  Res. 

Postfach  510153 

D-3000  Hannover  51 

FEDERAL  ItEPUBLIC  OF  GERMANY 

University  of  Hawaii 
Institute  of  Geophysics 
■iTTN :  Dr.  Daniel  'Walker 
Honolulu,  HI  9G1322 
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Processor  Peter  liarjes 
Inscicuce  tor  Ceuphysik 
Kltur  Universicy/Bochum 
F . 0 .  Box  102148  4630  Bochum  1 

FEDERAL  REPUBLIC  OF  CERMANY 

Professor  Brian  L-N.  Kenner c 
Research  School  of  Earch  Sciences 
Insticuce  of  Advanced  Scudies 
C.P.O.  Box  4 
Canberra  2601 
AUSTRALIA 

Or.  B.  Massinon 
SocieTe  Radiomana 
27,  Rue  Claude  Bernard 
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Dr.  Pierre  Mechler 
SuCieCe  Radiomana 
27 ,  Rue  Claude  Bernard 
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DARIA/ NMRU 

1-100  Wilson  Boulevard 
Arlington,  VA  22209-2308 

Dc  .  Pc  car  Basham 
Geological  Survey  of  Canada 
1  Observatory  Creseut 
Occowa,  Ontario 

Canada  kia  0Y3 

Dr.  Robert  BlanJiord 
DARPA/'NMRO 

1-',00  Wilson  Boulevard 
Arlington,  VA  22209-2308 

Saudi  a  National  Laboratory 
ATTN:  Dr.  11-  B.  Durham 

Albuquerque ,  NM  U7185 

L)r.  Jack  Evecuden 
L)  JGS-Lar  tliquake  Studies 
3-^5  Middle  field  Road 
Menlo  Park.  CA  99025 

U • S .  Geological  Survey 
a 1 V A .  Dr.  1-  dauks 
Nac'l  Earthquake  Reach  Center 
J-.5  Middle!  ie  1  d  Road 
Menlo  Park,  CA  1)4025 

Dr.  James  Hannon 
Lawrence  Livermore  Nac'l  Lab- 
P.u.  box  8u8 
Livermore,  C.l  94  550 

U.S.  Arms  Concrol  &  Disarm.  Agency 
ATTN:  Dick  Morrow 

Washington,  D.C-  20451 


IJaul  Johnson 

ESS-4,  Mail  Stop  J979 

L^s  Alamos  National  Laboratory 

Los  alamos  ,  NM  s7o4l> 


M*  •  Ann  K.c  i* C 
U.  wA LJ A /  .li’livU 

L'.ju  Wilson  houlovarLl 
Arlingtou,  VA  222U9-2308 

Dr.  Max  Roontz 
US  Ue^c  of  Energy/DP  331 
Eorrescal  Building 
iuoo  Independence  Ave . 

Wo sh i m* ton ,  D.C-  20585 
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Dr .  W  .  H . 

USGS- 

Office  of  Earthquakes,  Volcanoes, 
u  Engineering 
Branch  of  Seismology 
345  Middle field  Rd 
Menlo  Park,  CA  94025 


Dr.  William  Leith 

uses 

Mail  Stop  92B 
Reston,  VA  22092 


Dr  .  Richard  Lewis 
Dir  .  Earthquake  Engineering  and 
Ceophys-i.es 

U  .3  .  Army  Corps  of  Engineers 
Box  631 

Vicksburg,  MS  39160 


Dr.  Robert  Masse1 
box  25046,  Mail  Stop  967 
Denver  Federal  Center 
Denver,  Colorado  80225 


R .  Morrow 
ACDA/VI 
Room  5741 

320  21st  Street  N  .W . 
Washington,  D.C.  20451 


Dr.  Keith  K.  Nakanishi 

Lawrence  Livermore  National  Laboratory 

P  .0  .  Box  806,  L-205 

Livermore,  CA  94550  (2  copies) 


Dr  .  Carl  Newton 

Los  Alamos  National  Lab  . 

P  .0  .  Box  1663 

Mail  Stop  C335,  Group  E553 
Los  Alamos,  MM  67545 


Dr  .  Kenneth  H .  Olsen 
Los  Alamos  Scientific  Lab. 
Post  Office  Box  1663 
Los  Alamos,  NM  87545 


iloward  J  .  Patton 

Lawrence  Livermore  National  Laboratory 
P .0  .  Box  808,  L-205 
Livermore,  CA  04550 


Mr  .  Chris  Paine 

Office  of  Senator  Kennedy 

SR  315 

United  States  Senate 
Washington,  D.C.  20510 
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AFGSR/NP 

ivTTM :  Colonel  Jerry  J.  Perriao 

Bldg  410 

Bolling  AFB,  Wash  D.C.  20332-6448 
HQ  AFTAC/TT 

Attn:  Dr.  Frank  F.  Pilotte 
Patrick  AFB,  Florida  32925-6001 


Mr  .  Jack  Rachlin 
USC5  -  Geology,  Rm  3  C136 
Mail  Stop  928  National  Center 
Rea  con,  VA  22092 

Robert  Reinke 
AFWL/NTESG 

Kir eland  AFB,  NM  87117-6008 
HQ  AFTAC/TCR 

Attn:  Dr.  George  H.  Rothe 
Patrick  nFB,  Florida  32925-6001 

Donald  L  .  Springer 

Lawrence  Livermore  National  Laboratory 
P  .0  .  Box  308,  L-205 
Livermore,  CA  94550 
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Dr  .  Lawrence  Turnbull 
0SUR/1JED 

Central  Intelligence  Agency 
CIA,  Room  5C48 
Washington,  D.C.  20505 

Dr.  Thomas  Weaver 

Lus  Alamos  Scientific  Laboratory 
Los  Aimos,  Nl-l  97S44 

AFGL/SULL 
Research  Library 

Hanscom  AFB,  MA  01731-5000  (2  copies) 

Secretary  of  the  Air  Force  (SAFRD) 
Washington ,  DC  20330 
Orr'ice  of  the  Secretary  Defense 
DDR  a  E 

Washington ,  DC  20330 
HQ  DNA 

ATTN:  Technical  Library 

Washington,  DC  20305 

Director,  Technical  Information 
DARPA 

1400  Wilson  Blvd . 

Arlington,  VA  22209 

AFGL/XO 

Hanscom  AFB,  MA  01731-5000 
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AFGL/  LW 

Hanscom  AfB ,  MA  01731-5000 
DARPA/PM 

1400  Wi 1  son  Boulevard 
Arlington,  VA  22209 

Defense  Technical 
Information  Center 
Cameron  Station 
Alexandria,  VA  22314 
(5  copies) 

Defense  Intelligence  Agency 
Directorate  for  Scientific  & 
Technical  Intelligence 
Washington,  D.C.  20301 

Defense  Nuclear  Agency/SPSS 
ATTN:  Dr.  Michael  Shore 
6801  Telegraph  Road 
Alexandria,  VA  22310 

AFTAC/CA  (STINFO) 

Patrick  AFB,  FL  32925-6001 

Dr.  Gregory  van  der  Vink 
Congress  of  the  United  States 
Office  of  Technology  Assessment 
Washington,  D.C.  20510 

Mr.  Alfred  Lieberman 

ACDA/VI-QA' State  Department  Building 

Room  5726 

320  -  2 1  St  Street,  NW 
Washington,  D.C.  20451 


